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We study two alternative definitions of localized states in the lowest Landau level (LLL) on a 
torus. One definition is to construct localized states, as projection of the coordinate delta function 
onto the LLL. Another definition, proposed by Haldane, is to consider the set of functions which 
have all their zeros at a single point. Since a LLL wave function on a torus, supporting magnetic 
flux quanta, is uniquely deflned by the position of its A^^ zeros, this defines a set of functions that 
are expected to be localized around the point maximally far away form the zeros. These two families 
of localized states have many properties in common with the coherent states on the plane and on 
the sphere, viz. a resolution of unity and a self-reproducing kernel. However, we show that only the 
projected delta function is maximally localized. Additionally, we show how to project onto the LLL, 
functions that contain holomorphic derivatives and/or anti-holomorphic polynomials, and apply our 
methods in the description of hierarchical quantum Hall liquids 
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I. INTRODUCTION 



Fractional quantum Hall (QH) states are formed when a very clean two dimensional electron gas at low temperature 
is subjected to a strong magnetic field A defining feature of these liquid phases is a quantized value of the Hall 
conductance = ve^ j'K, where v is the filling fraction, i.e. the density of electrons in units of the density of a filled 
Landau level. The topological nature of these states is also reflected in fractionally charged quasi particles^, obeying 
fractional statistic^, protected edge modes, and a characteristic ground state degeneracy on higher genus manifolds^. 

At a semi-classical level, an electron in a strong magnetic field can be viewed as a localized distribution of charge 
with the size of the magnetic length Ib — y/hcfeB, and an orbital angular momentum / associated with the cyclotron 
motion around the guiding center position R. For v < 1 sl\ the low energy states reside in the lowest Landau 
level (LLL), where the maximally localized states have zero angular momentum with respect to the guiding center. 
Technically, the wave functions of these localized electrons are coherent states with a Gaussian charge density profile. 
It is rather natural to try to use these coherent states as a basis when trying to formulate a theory of the fractional 
QH effect; the cooperative ring exchange theory using coherent state functional integrals is an early example'^. 

With the success of the composite fermion wave functions describing the Jain series for v = p/{2pm ± 1), and the 
idea of hierarchical states built from successive condensation of anyonic quasi particles, a significant effort has been 
put into finding ansatz wave functions directly in the position basis^ ^. The coherent states have rarely featured in 
such discussions, although their usefulness for improving wave functions was stressed in an early paper by Girvin and 
JacP. 

The construction of variational wave functions was put on a more solid theoretical ground by Moore and Read^^l 
and Wen'^ who proposed that representative QH wave functions for a large family of states could be constructed from 
conformal blocks in a conformal field theory (CFT). This proposal was based on the deep connection between the 
Chern-Simons theories, which give a low energy description of the QH liquids, and certain two-dimensional CFTd^. 
The strength of this approach, compared with i.e. the composite fermion description, is that the topological properties 
are manifest, given cer tain w ell supported technical assumptions^'^. 

In two recent paper Ji^lIIl it jg argued that the natural interpretation of the conformal block wave functions, is as 
wave functions in a coherent state basis rather than in coordinate basis. In the simplest cases the two interpretations 
are equivalent, but the use of coherent states allows for a comprehensive understanding of the full QH hierarchy, 
including not only the prominently observed states in the Jain series, but also more exotic ones corresponding to 
alternating condensations of quasi particles and quasi holes. 

Interesting topological properties are manifest in the QH wave functions on higher genus manifolds, such as the 
torus, which has genus g = 1. Although several examples of QH wave functions on the torus are known, including 
the Laughlin states and the non-abelian Moore- Read state, there is no understanding of the hierarchy at the same 
level as on the pland^, or on the spherd^. An important step towards this goal is to construct and analyze coherent 
states on the torus, and this is the subject of this paper. 

In 1963 Glaubei^^ pointed out that there is more than one way of defining a coherent state for the harmonic 
oscillator. The most common definition is as an eigenstate of the annihilation operator a. Equivalently, a coherent 
state can be defined as a state with minimal dispersion in phase space. In the LLL, the analogy of the latter definition, 
is a state that minimizes the spatial dispersion of the guiding center coordinates of the Landau orbit. Alternatively, 
this can be viewed as the projection of a position eigenstate onto the LLlI^. 

The structure of the paper is as follows: On a torus with fixed boundary conditions, the continuous group of 
magnetic translations that is present on the plane, is broken to a group of finite translations. This group defines a 
lattice, with number of points, where iV^. is the number of states in the LLL. In his original paper on QH wave 
functions on the torus, Haldane proposed a set of coherent states that is naturally defined on this lattice^^. In Section 



III we study these states, which we will refer to as lattice coherent states (LCS). In Section pVJ we construct a second 



set of coherent states in direct analogy to those discussed in the previous paragraph; i.e. by projecting position 
eigenstates, which on the torus are periodic delta functions, onto the LLL. Since these coherent states are labeled by 
the guiding center coordinate, we will refer to them as continuous coherent states (CCS). In Section |V] we explore the 
localization properties of both the LCS and the CCS for tori described by an arbitrary modular parameter, r. 

An important issue in formulating hierarchical QH wave functions on the torus, is how to incorporate holomorphic 
derivatives, and this issue is discussed in Section |VT] There we find that a holomorphic derivative projected on the LLL 
becomes a sum of translation operators. We also discuss alternative ways of mapping higher Landau level functions 
onto the LLL. 

Finally in Section |VII| we give an explicit example of how anti-holomorphic components can be treated, by inter- 
preting these functions in the basis of CCS. We will in this Section find that the required ground state multiplicity 
is not manifest from the CFT construction, but that it can be extracted by projection onto well defined momentum 
states. 

For a review of coherent states see Ref. |5n]and|lT]and the comprehensive book by Perelomov^^ for details. Recently 
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Kovalski et al. has written several articles about different coherent stateP^'^Sl on the circle, sphere and plane. 

In the following section we will begin by a short summary of some properties of charged particles in a magnetic 
field. 

II. CHARGED PARTICLES IN A MAGNETIC FIELD ON A TORUS 

To fix the notation, we here summarize elementary facts about charged particles on a torus, and in a transverse 
magnetic field. In particular, we give expressions for the LLL single-particle wave functions. 

A. Magnetic translation operators 

We use the Landau gauge A — Byk, with the magnetic field pointing in the z-direction, B = Bi. We use 
units where both the cyclotron frequency and the magnetic length are set to unity, such that uj^ ~ eB/m ~ 1 and 
£b — \Jhl eB — \. A general magnetic translation by a displacement L, can be written in the Landau gauge as 
t(L) = exp[L • (V + lyx) — • (L x r)]. We use the complex notation t[oL + = t(ax + /3y), and consider a torus 
spanned by the lattice vectors Li = Lj;X and L2 = Lax + T^yY ■ Interpreting the torus as a parallelogram with the 
opposite sides identified, L^-, Ly are the width and height whereas La is the skewness. For translations along the 
cycles of the torus we have 

t{a) = e"^- 
t(i/3) = e*^^e'^^« 

t(a + i/3) =e*5"'3i(i^)t(a). (1) 

Hence, magnetic translations around the cycles of the torus commute only if L^Ly = 2'kNs, where Ns is an integer. 
The torus can thus be parametrized by two parameters: The complex modular parameter 

T= -^(La+*L,), (2) 

and the number of magnetic flux quanta Ng . A general boundary condition on a wave function "0 is given by 

t{TLM-A = e'1'^4,{z), (3) 

where the phase angles 0i have the physical interpretation of fluxes threading the two cycles of the torus. We will 
use complex coordinates z — x + ly although all wave functions will have a non-holomorphic Gaussian part. The 
boundary conditions are not invariant under magnetic translations since 

t{L^)t{a + iP)i^(z) = e'^'I'^+P^'^hia + iP)tp{z) 

t{TL^)t{a + iP)tlj(z) = e^-^^+zSLA-aLj^^-^ ^ il3)^(z), 

which follows from the magnetic translation operator algebra. Translating in one of the principal directions will 
change the boundary conditions along the conjugate principal axis. From this we see that the boundary conditions 
are invariant under a subset of magnetic translations, F = a + i/3 = + = j^n + jf^mr with integers n and 

771. It is convenient to use the notation 

Lx 27r 

Vn = —^n = — n (4) 

to parametrize the natural sub-lattice formed by these translations which preserve the boundary conditions. Note 
that XnUm = Stt^j^, so that e'^"^"^ = 1. We will call the displacement operators that move one A^^ith step in the 
principal direction 

h = t{xi) 

t2 = t{TXi) (5) 

respectively. In the following we shall fix the boundary conditions to 0i = 02 = 0, but all results can trivially be 
extended to arbitrary <j)i using the magnetic translation operators. 
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B. Bcisis states 

The Hamiltonian in the Landau gauge, with units restored, is 

1 2 1 / ,n 

h=^pI + i^{p. vf. 

Zni Am c 

A simple and instructive way to construct the torus wave functions is to first consider a cyHnder with the axis along 
the y-direction. The simultaneous eigenfunctions of ti and H , are easily obtained as 

X„..(z) - —L=e-^y'^H,Xy - 2/.)e-5(^-^=)' (6) 
V -t^^V"' 

where i?„ is a Hermite polynomial, and Hxn.s — tiujc{n + ^)xn,s- The eigenvalue of ti is t\xn,s = e~'^'^=x„,s. The 
states Xns can be obtained using ladder operators in the form = -^{y—ys~dy). Using the form of ^ one finds that 
dxXn,s = ~WsXn,s- We can build ladder operators that are s-independent as a) = -^{y — i2dz) and a — -^{y — i2dz)- 

These operators are also independent of the choice of (f)i, which a\ is not. The operators and a are also the ladder 
operators on the torus. Using a, a) and t2 we can, after fixing the fluxes 0^, generate the full set of torus basis wave 
functions. 

We will however here make an explicit construction using Xn,s and the required properties under tj. We want to 
build a linear combination rjn^s = 'Ylm r^"i,rXm,r with the properties: Hrjn.s ~ {n + ^)rin,s, t\'nn,s — e^'^'^°?7„.s and 

t2Vn,s — Vn,s-i- We must then choose am,r — Sm^nSi^"^ exp{t^yrUir)- Here Si^'^ = J2t ^r+N,t,s is a periodic Kronecker 
delta, and = jf^r is the skewness of the Xr x rxr lattice. In the LLL, the basis functions can then be written as 



t 



or in a more compact form as 



NsT . (8) 



Here 1^3 is the third quasi-periodic Jacobi theta function. 



k 

These functions are orthonormal on any translation of the fundamental domain, x L^t. In an orthogonal basis 
r]s diagonal in ti, i.e. tirjs — e^^^y^rjs, the operator t2 will act as t2?7s = ^7^+1 thus generating the full set of basis 
states. Since ti and t2 do not commute we may instead choose to diagonalize the t2 operator and would get the 
eigenfunctions: 



T 



(9) 



The t2 eigenfunction can also be obtained from r/s by performing the modular transformation r — > — - followed by 

Irl . ■ 

letting z — > ■'^^z and applying the appropriate gauge transformation. 



III. LATTICE COHERENT STATES 



All the states in the LLL are uniquely defined by the positions of the Ng zeros in the wave function, so we should be 
able to engineer a spatially localized state by choosing the position of these zeros appropriately. Haldane and RezayP^ 
proposed a candidate for such a localized wave function, obtained by putting all zeros in the same point. Because the 
continuous translations of the plane have been broken down to a discrete set of translations, generated on the torus 
by ti and <2i this point cannot be chosen arbitrarily. Instead we will get a family of wave functions. These wave 
functions form an over-complete set, and we shall refer to them as lattice coherent states (LCS). One usually thinks 
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of coherent states as states that minimize the spatial dispersion. As we will see in Section |V] the LCS are. in this 
sense, only coherent on a rectangular (5R(t) = 0) torus. These states are interesting since they naturally incorporate 
the finite translation structure of the translation operators that preserve the boundary conditions, ti and t2- To our 
knowledge this is the first time these states are carefully examined, so we shall work out their properties, including 
in particular, formulas for the resolution of unity, and the reproducing kernel, in some detail. 
The most general wave function in the LLL on a torus is 



11 i-iz-^j) 



(10) 



where is the position of the j:th zero. Here di is the first Jacobi theta function. 



di{z\ t) = 2 ^(-l)"e*^^("+5)' sin ((2n + 1) z) . 



By demanding that ^'(-z) obeys periodic boundary conditions defined by (|3|, we get relations on k and S, = ^ X]j 0- 
These can be written as 



ikL^ _ (^_2)"'^= 



(11) 



Solving (111 leads to the relations 
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^ Vr 
k 

Lr,, 



(La — Lx) 



-k. 



Here n and I are arbitrary integers. We can thus write as — xi[m ~ nr] 



[r — 1], where m = n + I. We 



define zj = £,j + -^{^ + t)Lx'; since we expect that the LCS, where all the zeros are at the same point, will have 
the maximum at the position diametrically opposed to ^. These coordinates, Zj, will give possible positions of the 
maximum value. Under this re-parametrization the general wave function can be written as 



2 1. We / 

Lx ^ 



(12) 



To proceed we will use the Jacobi theta function identities 



"■(-I 



±^2{z\t) 



to transform -dx into dj,. Here -Qi is the second Jacobi theta function 



§2{z\ ^) = 2 ^ e^^r{-^+\f cos ((2n + 1)) z) . 

n=0 



We have up to a scale factor and a constant phase 



^iz)^e-'^e^y'^^Y[dJ^{z-z,) 
j=i \ ^ 



(13) 



Defining the mean value of the zeros, z = '^j=i ■^i ^"^'^ assuming all Zj ~ z we get the requirement z ~ Xm — XnT, 
where again n and m are integers. Changing n — >■ —n we get the full wave function, with some n and m specific 
normalization J^nm j to be 



(14) 
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where 



iT. For numerical evaluation, equation (14) is a useful expression. However, for analytic 



manipulations this is not the best way of writing 4'nm{z), and it also leaves unanswered the question of how to 
calculate the normalization A/'nm . We solve for the relative normalization by using the magnetic translation operators 



to transform ipoo into 4'nm- We see by inspection that \Afn. 



where A/" = A/ni 



A. Fourier expanding ipnm 

To bring ipnm to a form that facilitates further manipulations we expand "d^" in Fourier modes as 



1 



oo 
K=~oo 

In the second line we have made the substitution kj — ^ + kj, where K = such that kj is the deviation 

from the mean value of kj. We hide the nontrivial sum over k in the value Zj^, 

oo 

Zk^ J2 e""^."=i^?. (15) 

The constant Zk can, together with the factor e ^'^'^ , be interpreted as the partition function of Ns particles on a 
circle with the total angular momentum K. The expression for "000 is thus 



0oo(z)=A/- J2 Z;^e-5(2/+!^^)'e*^'^("+^) 

which for general ipnm becomes 

oo 

^nm{z)^Af ^K+„e-3(^+^^)%'^^("-=^'")e'5''^"^. (16) 

We note that for equations (16) and (14) to match, we must have Afnm = e*''"^^'"+5^"'^^7V. The price for this simple 



formula is that we now have a set of complicated constants Zk- Fortunately for our purposes it is sufficient to know 
the periodicity property Zk+n, = Zk- It will turn out that the details about Zk can be hidden in the normalization. 
Zk is periodic because it only contains deviations from the sum mean value 



B. Overlap and normalization 



The overlap between the two states ipnm and ipn'm' can be calculated by noting that the a;-integration gives Kronecker 
deltas, that allow us to combine the y-integral from an incomplete to a complete Gaussian integral. We decompose 
the sum over K as K = I + tNg, such that X]k=-oo ~ St^-oo- '^-'^^ overlap on any fundamental region is thus 

Ns oo „L„ 

;=i t=-oo-^o 
After the Gaussian integral is performed the overlap reduces to 

ii^n'm' li^nm ) = V^L^M^ ^ Zz+„Z,Ve^^' . (17) 

1=1 
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Choosing m! = m and n' = n we get 



=1 



which defines the normalization constant in equation (16 1. 



C. Resolution of unity 

Just as the usual coherent states on the plane, the states ij^nm allow for a simple resolution of unity in the LLL 
given as 



VlLL = ^ X! 1^"'") (V'nml • (19) 



m,n— 1 



We now prove this relation, which will be useful in applications of the coherent states. First notice that from ( 16 1 we 
have 

tl ^2 I'/'nm) — e ^" \lpn—n',m—m') (20) 
(V-nml tft^' = e'«"'(-'"+-™') {i,„+„>.rn+m'\. (21) 

and from this, it follows that 

by changing the summation index. Since ^^^^ ^ IV'nm) (V'nml commutes with any translation it must be proportional 
to the identity, i.e. ^ IV'nm) (V'mnl = cT'lll where c is a constant. To determine this constant, it is sufficient 
to calculate a single matrix element. To do this we first express Vlll in the complete basis \ris) given by the wave 
functions ([8|, such that 

(z = IT'lllI z = 0) = ^ |(z = |2 = -^-i— e~^^'-^y'+^'^y^'^e~^^y''-'^«^'-^y' . 

s r,t 

Calculating the same matrix element using the lattice coherent states, we obtain 

^ |(^ = |V„™) I' = iV, • (z = \Vlll\z = 0) . 



This implies that c — Ng and thus proves (19 1 



D. Self-reproducing kernel 

Since the, un-normalized, coherent states fwiz) on the plane can be obtained by projecting a position eigenstate 
onto the LLL, they are closely related to the holomorphic delta function. The state (^u, is in fact a self-reproducing 
kernel on the space of LLL wave functions, i.e. J (Pw Lp^^{z)ip{w) — ip{z). If we instead convolute the kernel ip^iz) 
with a wave function that has components in higher Landau levels, these will simply be projected out. It is not 
obvious that there is a similar kernel for the LCS ipnm, but we can in fact prove that 

^P{z)^S~^ e'^?/"""^,„„(z)V;(z„„0, (22) 

m,n— 1 

where ipi^z) is again an arbitrary wave function in the LLL and Znm = + XnT. If we substitute 'ijj{z) with an 



arbitrary wave function (j){z, z*) then (22) defines a map, PlcSj from the space of all wave functions to the LLL wave 
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functions. Although this map is not a projection on the LLL. it might still be useful in constructing quantum Hall 



wave functions, and we shall comment on this in Section VI 



We now prove (22) by first looking at the specific case where i^{z) is a basis wave function. With a bit of algebra 
we can establish 



m,n— 1 

where rjs is any basis wave function defined in From this follows 

Ns 

^(Z,fc)=5-1 ^''^''"""i^nmizMZnm), (23) 



which holds for all l,k £ Z where is any state in the LLL, and 5* = Ng^QQ{0) is a normalization factor. That (23 1 



holds for Ns x Ns points is more than enough to ensure that (22) holds for arbitrary z. To prove this, consider the 
difference between the left hand side and right hand side of (23), for arbitrary z. This difference, ^(z), must have 
zeros at all z = Znm- But ^(z) can fulfill the boundary conditions, if and only if it has the form of ( |10[ ), which only 
has A'^ zeros. Therefore ^(z) — identically, which concludes the proof. 

This LCS-map should be used with care, because it is not a projection onto the LLL. Using T'lcs on non-LLL 
function will in general give nonzero results. This is seen by considering 5{z — z') which has components in all Landau 
levels. It is obvious that Vi^csS^z — z') will give zero even though we know that S{z — z') has components in the 
LLL. Thus the effect of T'lcs is that the contributions from non-LLL states precisely cancels the LLL part, except at 
z' = Znm for which the contribution is divergent. 

It is instructive to consider Plcs for the special case of Ns — 1. In this case there is only one state and Plcs^I-z) = 
^°°|q| 0(O). If (/) is a basis state ??n,o(-2) we can, via a simple parity argument, show that ?72n+i,o(0) = and r?2n,o(0) ^ 
in general. For iV^ 7^ 1 a similar analysis is more difficult, since r\n^s should be evaluated analytically in iV^ x Ng 
points together with tjjnmizi^p). Numerical studies suggest that 7^LCS^/2n+i,s = and 7'LCS^2n,s 7^ does hold for 
general Ng. This would mean that if (p is restricted to the two lowest Landau levels, then Vi,cs4> = 'Plll4'- 



IV. CONTINUOUS COHERENT STATES 



In the previous section we introduced the LCS wave functions as candidates for coherent states. From the point 
of view of describing localized states, the LCS have the drawback of only defining states that are localized around 
the finite set of points spanned by the lattice vectors xi and txi. It is complicated to construct a localized particle 
around some other point with the LCS. This brings us to the notion of continuous coherent states (CCS) that are 
obtained as the projection of a position eigenstate on the lowest Landau level. We thus define 



Vwiz) = 7'lll'5(z - w). 



(24) 



Here w = 

Section [III[ The overlap between two coherent states ip^^ and ip^^ is readily obtained as 



ly' . As a projector we can either use some basis states Plll = X^s I'?'*) or we can use (19) from 



{u\w) = 2j'/*(^)'7s('^) = Vwiu) 

S 

From the definition of (pw{z) also follows a resolution of unity 



(25) 



ip{z) ^ d wipyj{z)tp{w) 



for states in LLL and zero otherwise. For analytical manipulations we can either work with (19), and construct the 
LLL CCS as 

^wiz) = ■^Ylpnni{z)^p*„,{w), 
^ mn 

or work with ris{z) and obtain ipw{z) using (pw{z) = X^fl^i Vs{'^)''ls{z)- We can note that when we work with arbitrary 
boundary conditions, given by (|3|, we can use the transform 
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where 7 = §^((^2 + ^it) and t^^^ acts on the z coordinate. After some algebra, where we use the sums over m and n 
to cancel the Zk factors against the normalization, we get 

^ ^ K,t 

With the aim of getting an expression in terms of "iJ-functions, we rewrite tp^ as, 

ipU^)^P{y,y')xYiz,w*), (27) 

which factorizes into a Gaussian part 



and a holomorphic part 



where Et can be simplified as 



Y{z, w*) = e-™**-^''e~'''^=^*'£;t, (28) 



Et^^3[ ^{W* - z)+TTtf 



^ *3(r) ) . (29) 



N. 



In order to have an expression for Y(^z, w*) that is compact, we would like to find a transformation that removes 
the term irtf from the i^s-function. However since 1^3 has only 3(t) as second argument this is not possible. In the 
rectangular case where ^(t) = we can proceed. 



A. Purely imaginary r 



In the following we will assume that t is purely imaginary i.e. t = ijf-. Even so, we will obtain different functional 
forms depending on whether is even or odd. This is related to the structure of the zeros. For an even Ng the zeros 
will be divided into two groups. Each group of zeros will, for any w, be regularly spaced on some line parallel to one 
of the fundamental axes. For an odd Ng the zeros cannot be divided into two grops, and the distributions of the zeros 
is more intricate. For even Ng we will have (31 1 and for odd Ng we will have (32). It will be convenient to define: 



T- 



L. 



-{z~w*) 



T+ = —{z + w*). 



(30) 



In the case of an even Ng we can use the i^-function relation 



^■i{z + krirl r) — e 



-2ikz 



which reduces ( 28 1 to 



Y{z,w*)^d3[ T 



2t\ / Ng , 
Ng n 2 



NgT 

2 



The full form of the continuous coherent state for even iV, is thus 



e-W+y'^) 



m T 



Ng n 2 



NgT 

2 



(31) 



We can see that the zeros of this coherent state lie on two lines intersecting sX z = w + }^[\ + t)Lx as drawn in the 
left panel of Figure [T] This structure is entirely different from the LCS, where all the zeros are at the same point, as 
indicated in the right panel of Figure [T] 
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CCS structure Df zeros LC5 structure of zeros 




u u 



Figure 1. Structure of zeros of continuous coherent states ip^ (left) and lattice coherent states -^nm (right). The red circles 
represent the maximum density, whereas the black circles represent zeros. On the rectangular (5R(r) = 0) torus for Ns even, 
the zeros of the CCS form a "cross" centered over z = w + + t). For general r and/or odd Ns the structure of the zeros 
is not as easy to describe. 



For odd Ns, the expansion of Et in (29 1 depends on whether t is even or odd. If t is even we basically have the 

2 



same result as in the case where N^ is even since we know that is an integer. For the odd i-values we will need 



to use some more ?9-function relations 



^^3 ( 2 ± ^7I"T 



T =q--^e^"d2{z\T) 



d2{z + m:T\T)^q-'' e~^"'^^2iz\ t) 
such that is transformed into i?2- The final result for ipw{z) in the case where Ns is odd is 



2t 



dj{T+Ns\2NsT) 



(32) 



B. Generic r results 



In the generic case we see that we cannot simplify ifyj to a finite sum of t?-functions. We get for even Ns, after 
shifting away the 7rt3(f) part, 



e 2 tf 2 a 



f"^(")*'i?3( T- -7rt3fi(r) 



N. 



i3(r) 



(33) 



We see that if it was not for the term 7r<5R(T), in the i?-function, the projection would be given by (31 1 albeit with 
T — i5(t). For odd Ns we get a similar result resembling (32|2Si. 



V. LOCALIZATION BEHAVIOR OF CCS AND LCS 



In this section we will analyze the localization properties of the LCS and the CCS wave functions by calculating 
OxCTy for different t and Ns, where cr^ = {x^) — (x)^. We will for simplicity consider projections of w = unless stated 
otherwise. We will calculate the spatial dispersion using {A{x,y)) = J dx dy A{x , y) \(pw(z)\'^, where the integral runs 
over a fundamental domain il, chosen so that (x) and (y) are at the center of fl. For highly skew tori the LCS might 

not be localized around only a single point. For instance, we will see that if we set r = «^ + |, which corresponds 
to a triangular lattice, then ipnmiz) develops two maxima. These will by symmetry be located at the center of each 
triangle and form a honeycomb lattice, as can be seen in Figure |2] This feat ure should be stable even for Ns oo, 
meaning that the LCS can never describe a localized state if r is tilted enough. 

Figure [2] shows how the spatial profile changes as r is tuned away from 3?(t) = 0. The CCS retain their circular 
shapes even at small Ns- For the same change in r, the LCS become very distorted^^'. The LCS distortion is a 
combination of a geometric effect, and that the LCS develop multiple maxima. The geometric effect is because fir 
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Figure 2. The spatial profile of CCS (upper) and LCS (lower) at Ns — 4. No axes are drawn as the black square shows the 
fundamental domain, fi. The columns correspond to r = + Tx where a) = Q, h) Tx — \, c) Tx = \ and A) Tx = \- The 
values of Tx correspond to a rectangular, two general and one triangular lattice. Lighter colours correspond to larger values 
of I'i^l'^. The CCS nicely reshapes itself whereas the LCS becomes very distorted and develops a double maxima. The double 
maxima should be stable even for large values of Ns because of symmetry. 



(TjfO-y for Ns — 3 and + r y 




a.Q 0,2 0.4 06 o.a 1.0 1.2 iA 5 10 15 20 25 30 35 

a b 

Figure 3. The spatial dispersion for r = i measured as Gx<Jy. Left: OxOy for with variation of w. The smallest dispersion is 
at = whereas the largest is at w = |(1 + r). Right: Ox(Jy for a CCS, (red) on w = Q and (green) on w — |(1 + r), and LCS 
(blue) for some values of A^s. The CCS in general have smaller dispersion than the LCS but the dispersion of CCS depends on 
w. For Ns = 1, 2, 3 the value of axCJy coincides for LCS and CCS because the wave functions are the same these cases. 

has its center at the boundary of fir+i, if they share one corner. Because of this effect we need to be careful with the 
measure ax<Jy It is not guaranteed that the maximum of is located aX z — (x) + i{y). 

Figure [3] shows how the dispersion depends on the number of fluxes, Ns- The dispersion is signiflcantly smaller 
for small values of Ns < 10, because at these sizes the torus is so small that the CCS have the same width as the 
torus. At Ng > 10 the CCS have saturated at the dispersion value expected on the plane, whereas it takes the LCS to 
Ns > 40 to reach the same values. We can get a feeling for why the CCS are more localized than the LCS by noticing 
the structure of zeros for the CCS and LCS. It is likely that as the LCS have all their zeros diametrically opposed 
to the maximum value they do not localize the wave function as efficient as the CCS do where the zeros are spread 
along the border of the torus. 

It is noteworthy that the axCy-vahie of the CCS depends on the precise value of w. If w — x„ + rxm the dispersion 
is at a minimum, but if w — x^^^i + tx^^i the dispersion is at a maximum. This suggest that the xi x xit lattice 
is encoded in the CCS. In the left panel of Figure [3] (T^ay is plotted for r = i and iV^ = 3 where w is varied over 
the fundamental domain. A version of this is also seen in Figure H where in the upper panel the contours of \(pw\'^ is 
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Figure 4. The spatial profile of CCS at A^'s = 3 for = y f » where a) w = 0, b) «; = jXi and c) w = ^xi. These correspond 
to a rectangular lattice where we move the w away from w = 0. No axes are drawn as each figure is one fundamental domain, 
n. Notice how the spatial profile changes as w is tuned away from w = 0. This is a consequence of the zeros of Lp^ moving 
around. Upper: Contours of \'Pw\'^ with fl centered at 2 = w. Brighter colours correspond to larger values of |</Jiup. Lower: 
Contours of log \ (pw\ with centered at z — |(1 + t)xi. A logarithmic scale is used to make the positions of the zeros (filled 
red circles) easier to see. 



shown for w = 0, w 



'At — i and r 



4«- 



|a;i and w 

For small values of Ng the dispersion for w ^ Q can be higher than that of the corresponding LCS. However, already 
at A^s = 9 the maximum and minimum dispersions for CCS at different w are practically indistinguishable. What 
happens at a more technical level is that the zeros of the wave function start moving around as can be seen in the 
lower panel of Figure |4] Here the torus region is held fixed, centered at r = ^(1 + t)xi, to make it easier to see that 



the sum of all the zeros is fixed as required by the boundary conditions in equation (111 



Figure [5] shows what happens when r is tuned away from the square lattice st t — i. We see here that there is a 
wide region of r where the CCS have smaller (JxCfy than the LCS. However at large values of 3(t), corresponding to 
Ly > L^, the LCS seem to have a better dispersion than the CCS. One should mention that this effect is tiny and 
starts to show first when the torus width is so small that the CCS wave functions would be wider than the range. 
If we instead change the real part of r then the CCS functions become more localized in comparison with the LCS. 
What happens with the LCS is that as r — >■ t + ^ the torus geometry becomes triangular and the LCS develop two 
distinct maxima, as mentioned above. 



VI. HOLOMORPHIC DERIVATIVES PROJECTED IN THE LOWEST LANDAU LEVEL 



QH wave functions that describe hierarchical states on the plane and on the sphere, are most simply expressed 
in terms of holomorphic derivatives acting on holomorphic polynomials. In the composite fermion approach, the 
derivatives are remnants of the anti holomorphic coordinates z* present in higher LL wave functions^ . In the approach 
based on CFT, they appear since the electron operators typically are Virasoro descendants of primary fielda^^. On the 
plane the derivatives are in principle easy to handle, although it can be numerically difficult to take many derivatives 
of the polynomials, which are of order N'^, where N is the number of particles. On the sphere the situation is more 
complicated because of the finite dimensionality of the space of LLL wave functions, but recently the full hierarchy 
has been constructed also in this case—. On the torus the situation is more intricate. As discussed in Ref. [3S1 in the 
context of the CFT construction, the electrons are still most naturally represented by descendants of primary fields. 
The corresponding conformal blocks, which are well defined on the torus, contain holomorphic derivatives. Although 
these blocks do not satisfy the proper QH boundary conditions, for the simple case of the Laughlin states, a subspace 
of functions that do can be found, resulting in precisely the torus Laughlin wave functions, constructed by Haldane 
and RezayP^. For the hierarchical states, which on the plane involve holomorphic derivatives, this procedure fails. A 
partial solution to this problem was proposed in Ref. 1281 where it was noticed that if the holomorphic derivative was 
replaced by a certain finite translation operator, the wave functions were well defined, and had good overlap with the 
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o'x^v for Ns = 10 and r=j + Tj^ ^x^v fo^ = 10 and r=/ Ty 




Figure 5. The dispersion a^o-y for CCS (red curves) and LCS (blue curves) at Ns = 10 as r is varied, a) Here r is varied 
from r = I to T = I + 1. As r changes the CCS stay equally localized whereas the LCS becomes strongly distorted, b) Here 
T is varied from r = i to r = 16z. As r is changed there is a region where the LCS have slightly better dispersion. For large 
imaginary values of t then axCy go to zero roughly as J: — , which is expected when Lx is smaller than £b and the dispersion 

in Lx direction is suppressed. 



coulomb ground state, at least in certain geometries. 

In this section we shall approach the problem from the point of view of coherent states, and show that the occurrence 
of finite translations in Ref. |28iis not a coincidence. The finite translations are related to how holomorphic derivatives 
of torus LLL wave functions, which are not themselves in the LLL, behave when projected back to the LLL. 



A. How to project holomorphic derivatives onto the LLL 

There are several routes to derive expressions for the projected derivatives. The easiest way is to use ladder operators 
and cylindrical wave functions and their overlaps as defined in Appendix [B] It is important to note already now that 
the result will be sensitive to constant shifts in the gauge potential. Constant shifts in the potential A — > A + y' , 
correspond to translations of the fundamental region. An easy way to see that the choice of fundamental region 
matters is by considering ^(z) and cjj^z + rLy), where is a LLL wave function. We see that dx4>{z + iLy) — 
—iLye~'^^^^ 4>{z) + e~'^^^vdx4>{z)- Since dx(j>{z + iLy) ^ dx4>{z) the following projection will not be the same. 

We start by examining the effect of acting with one derivative on the basis states in the LLL. Using the ladder 
operators we can express how the operators y, dy and dx act on 770, s = Vs- Since both y and dx break the torus 
boundary conditions, their actions are best described in terms of the cylinder functions Xn.s- Letting y, dy and dx 
act on the basis functions rins we get 



yvn 



dx^ns — ^^ns 

where A„s = '^^ys+tNsO'n.s+tNsXn,s+tNs- The coefficients a„^s were defined in Section hll Note that that A„s has 
components in all Landau levels, since Xn,s is not in a single Landau level. The holomorphic derivatives are obtained 
as 

dzVns ^ ^ (^Kis + ^^('^ ^ a))rins^ ■ (34) 

As long as we are in the LLL then (?7o,s l^z] 77o,r) — {vo,s\dz\'r]o,r)i since (770,3 770, r) = 0, which means that 
to investigate the effect of dz we must calculate (770, r \dz\ Vo,s) — Yi (^o,r |Ao,s)- For this purpose we need to know 
the overlap of the cylinder functions evaluated on the torus, and this calculation is performed in Appendix B. It 
is important that (r/o.s l^^zl %,r) is proportional to 5^ g with an s-dependent constant Gg. The effect of projecting a 
derivative down to the LLL is 

T'lllS^ |7/o,s) = ^Gs ho.s) , (35) 
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so that Gs+N^ = Gs and is thus periodic. This enables us to write 

where ai — e*^=^'Gs is the discrete Fourier transform of Gg- Because the result does not depend on s it holds 

for any state in the LLL. and consequently 

for any LLL state. 

The same result can of course be obtained using LCS or CCS, but the calculation is more complicated^^ 



B. Higher order derivatives and derivatives in the thermodynamic limit 



In (34 1 we saw that dzrjns will have components in all Landau levels because \ns- It is therefore difficult to give a 
closed formula for ViAA^d^rjns- We can however give a general formula for {ViAA^dzY (j){z) since at each projection we 



can apply ( 36 ) . The effect of multiple derivatives on one state can thus be expressed as products of sums of translation 



operators, which we can rewrite as a single sum of such operators: 



With the definition above, and ai''^ = ai, we get a recursive relation 



N. 

n " s 

where the coefficient am is simply the m:th Fourier coefficients of G^. We cannot calculate the coefficient ai ana- 



lytically, since it contains limited Gaussian integrals (see (C2| in Appendix O), but we can take the thermodynamic 



limit where both Ng and Ly become large. In this limit take almost continuous values and we can approximate 
G, with a continuous function 



G{0=^V^-^ E / dye-y- (38) 



k— — oci 



such that Gs = G{^)^. By construction G(^) has periodicity 1, where we specialize to a torus centered around 
z = 0. We can simplify the expression for G(^) by first taking a derivative, then expanding in Fourier modes and 
finally integrating in ^. The result is 



n=l 

and the Fourier coefficients of this function are 



G(0 = E sin(2^<). 



(-l)V 



(39) 

yvi/t 

The n = 1 term will dominate in the Lj, — > limit where G(^) becomes a sine function, which as Ly oo is changed 
to a sawtooth function. At first this gives the impression that the dominance of the first term is diminished in the 
thermodynamic limit. This is certainly true but of little consequence since the number of terms that are reasonably 
large scales only as niargc ^ The relative number of translations thus decreases as ^issE ^ go that in the 

thermodynamic limit the relative number of the relevant translations goes down. This is actually what one would 
expect since in the LLL, possition cannot be specified closer than one magnetic lengtlPfll 



15 



C. Holomorphic derivatives and LLL wave functions 



We now return to the question of how to treat the holomorphic derivatives in the LLL CFT wave functions. The 
problem we face is that from CFT we obtain wave functions that naturally contain derivative factors dz- Since these 
derivatives violate the boundary conditions and force the CFT wave function into higher Landau levels, we would like 
to be able to project the dz onto the LLL in a consistent manner. We have above found that under projection to the 



LLL a holomorphic derivative transforms as 9^ ^ I? = ^ ^ o.ii''xi where a; where given by (37 1. This suggests that the 
substitution 9^ — > I? in the relevant conformal blocks would yield good hierarchical QH wave functions on the torus. 
This was in fact proposed in Ref. [25] but without any motivation, and with the coefficients a; as free parameters. 
This simple prescription will however not work for the following reason. The quantum numbers for the many-particle 

wave functions can be taken as, Ti'ipj{zi, . . . , zpf^) = e ^"^ipj{zi, . . . , zn^) and T2tljj{zi, . . . , zn^) = ^j+i{zi, . . . , zn^), 
where Ti = Ylk *i,fc ^^'^ ^2 = Yik ^2,fc translates all coordinates one ti-step or i2-step respectively. It is straight-forward 
to see that T> does not commute with neither T2 nor Tj* , and consequently does not respect the quantum numbers. 
What will work, is the substitution, 

n5,^^(a,Tf + 6fcr2'=), (40) 

j k 

and it was in fact an ansatz of this form that, for reasons of simplicity, was tested numerically in Ref. .28, Since the 
simple substitution, dz ^ T) does not work, it would be desirable to have some other guiding principle for finding the 



coefficients and bk in the expression (40 I, and it turns out that modular invariance is very important 

We now turn to the remark concerning the lattice coherent states. In Section |III| we defined the map Plcs from 
the whole space of torus wave functions onto the LLL. This map could be used to consistently bring the CFT wave 
functions downs to the LLL, even where there are arbitrary functions of dz- The states that we have in mind are 
the ones that cannot be described in the hierarchy picture by simple condensations of quasi particles. In these states 
anti-holomorphic components naturally arise as i3i(zi — Zj \ r). Transforming z ^ dz ^ T) would be intractable given 
the amount of terms, but using T-'lcs we could simply evaluate the trial wave function on the multidimensional lattice 
{z-minn ■ ■ ■ , z^nj^ un } ^^id construct the LLL counterpart from there. 



The usage of Plcs also opens up for an application related to calculating overlaps. Given (22 1 we can calculate the 
overlap between two LLL states, i/j and (j) as 



(0|<^)=5 ^ 0*(z„/„.)((9(z„™) (m'n' |mn) , 



-2 



where we here restrict ourselves to 5ft(r) = for a cleaner notation. For an iVg-body state (zi, . . . zjv^), which can 
be expanded as 



n,m J^l 

the many-particle overlap is 



n',m' n,m j = l 



,t4JV, 



As this stands we would need to evaluate cj) and ip at approximately ^ , points each and calculate the ,^2 cross 
terms. If we are bold and assume that only the diagonal {nm \nm) terms will contribute, and all other will either be 
small or be averaged to zero, we can use the approximation (rt'm' \nm) w 5n.n'Sm,m' ■ In this case we would get 



S " ^ ^ 4' (^nimi J ■ ■ ■ I ^njVe'WJVe ) (^nimi 7 ■ ■ ■ 7 ^rijv^ mjv^ ) 



which gives us roughly ^ , points to evaluate. This is still a large number but we can here hopefully use that the QH 
system is strongly correlated, to remove the (majority) of terms that are almost zero because two electron coordinates 
will be close to each other. This could be used, not only as a systematic way of truncating an overlap calculation in 
the LLL, but also to lessen the computational burden in using T'lcs to map higher Landau level wave functions down 
to the lowest one. 
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VII. USING CCS TO PROJECT TO THE LOWEST LANDAU LEVEL 



In this section we give an example of how the torus coherent state kernel can be used to project many-body wave 
functions onto the lowest Landau level. Following Girvin and JaclP we consider a short distance modification of the 



1/ = - Laughlin wave function 

q o 

q+2p 

4q 



i<j i<j 

which has both z and z components but no external derivatives. Note that due to the Landau level projection there 
is no simple plasma analogy, but ^'i still has the same maximal angular momentum as the usual Laughlin state. 

Since 5* i is only a short distance modification of Laughlin's wave function, we expect that this state is in the same 
universality class as the Laughlin state, and a good trial wave function. In the CFT approach to the quantum Hall 
hierarchy, given in Refs. 14 and 15, the wave function in the coherent state basis is given by products of holomorphic 
and anti holomorphic blocks as in (41 1. The wave function in this basis is calculated as the correlator of the electron 
operator V{z) = e*^'?+P^i(^)+*VP'^2(z)^ ^^le electronic wave functions is obtained by projecting on the lowest Landau 
level, which amounts to a convolution with a coherent state kernel. On the torus, the latter point of view is fruitful 
since the coherent state wave function, which is comparatively simple, obeys the same boundary conditions as the 
much more complicated projected electron wave functions. This is apparent since t{a) commutes with T^lll which is 
a convolution with the coherent state kernel, / (f^ ipiO- From this it follows that the same periodic boundary 
conditions that apply for the coherent state wave function t{Lx)ip(z) = e'^'^ip{z), also apply for the LLL projection. 
I.e. t{L,)VLLL^{z) = e'^Vi^i^L^z). 

On the plane we can factor the electron operator into a chiral and an anti chiral sector and directly evaluate a 
correlator to get ( [4T| ). The situation on the torus is more complicated since there are several chiral and anti-chiral 
sectors, and there is no procedure to factor the operators to directly get the trial wave functions. Instead we follow 
Ref. 28 and use the transformation properties of the different conformal blocks to arrive at an appropriate subset 



of the full correlator. The torus counterpart of the wave function (41 1 can be written as a single correlator of the 
full electron operator V{z) = e'v'9+P¥'i(z,2)+«\/P¥>2(z,z)^ where (pi, (p2 are compactified massless scalar fields with radii 
Ri = -jS= and i?2 = ^. 

For notation, and the technical procedure for constructing hierarchy wave functions on the torus, we refer to Ref. [28] 
and only quote the result here. The full correlator can be written as a sum over conformal blocks (Oiii ^(^i)C'b 



NiT)J2EuE2 ^Er,E2i^E^,E2 ^here 

^Ei.E2 = 



<;+2p ■ 

2q 



\q+p 



M- 



ffTE,.E2 {Z\t) 



^^y^\{d,{z 

i<j 

and N{t) is a normalization that will be of no importance here. Here Z = ^Ylik 



(42) 



-(zi — Zj). The sum 



over E\ and i?2 runs over a lattice spanned by Ej 
The conformal blocks are given by 



Ri 



and Ei 



m j R 



R, 



J'Ei Eo {Z\t) = e"VEl-fEl]^2^t[E^./^Z-E2yrpZ] 



- where Cj and rrij are integers. 



(43) 



Fixing the boundary conditions in the and tLx directions amounts to selecting a subset of the Ei, E2 lattice. 
This subset is parametrized by {Ei,E2) ~ k{y/q + p, y/p) + Tq where fc € Z. The simplest choice of the offset 
Lq — + P, y/p) gives the center of mass function 



{r+k)y'q+p,(r+k)^ 



{{q + p) Z - pZ\t {q + p) - Tp) 



(44) 



Here 



{z\t) is a generalized function such that Hr+i = Hr- Under center of mass translations T2 = Ylk ^ 



k ''2,fe 



we find that Tir — >■ ^r+i ^^'^ have recovered the expected 5- fold degeneracy of a = ^ state. The center of mass 



piece (44 1 implies one-particle periodic boundary conditions with phases 0i — T:{Ne — 1) + 2Trqr and (/)2 = 7r(7Ve — 1). 
We can select arbitrary boundary conditions (pi and 02 by translating the center of mass coordinate Z appropriately. 
After such a magnetic translation the final modified Laughlin wave function for (f>i — (j)2 — is given by 



q+2p 

2q 



\q+p 



a 



{{q + p) Z - pZ\t {q + p) - rp) 



(45) 
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where n = 1, . . . ,q enumerates the diflFerent momentum sectors and a = \{Nf. — 1). Furthermore, since we already 
have showed that the boundary conditions are preserved under projection onto the LLL. we have constructed a torus 
version of the wave function (41 1. 



We end this section by numerically comparing (45 1 with exact diagonalization of the Coulomb potential at r = i. In 
their original work Girvin and Jach noted that the Laughlin state ijM'^^ could be improved by considering components 
with p 7^ 0, and on the torus we observe the same thing. For p ^ Q then tJj^'^'P') is not entirely in the LLL. but 
the projected wave functions ViAA.ip'^'^'^^ do still have good overlaps with the Coulomb ground state. For A^'e = 3 
electrons and = 3 x 10^ Monte Carlo points, the {q,p) — (3,2) state has the best overlap with exact Coulomb, 
I (V'Couiomb IV'^^'^^ ) ? = 0.9999(4 ± 6^221. This should be compared to Laughlin, which has | (^Acouiomb l-iALaughiin) ? = 
0.9990(0 ± 2). For A^e = 4 electrons and A^ = 3 x lO'^ Monte Carlo points, the {q,p) = (3, 1) state matches Coulomb 
best, with I (V'Couiomb |V'^^'^^ ) ? = 0.9976(5±6) compared to | (V'Couiomb IV'Laughiin) P = 0.9792(8±3) for the Laughlin 
state. Going to larger system is exponentially difficult as the projection on the LLL requires an overlap calculation 
with all basis states. 



VIII. SUMMARY 



In this paper we studied two alternative ways of realizing coherent states in the lowest Landau level on a torus. We 
explored the set of lattice coherent state wave functions and found that they can be used to define a map Plcs which 
for the two lowest Landau levels acts like the true projection Plll- 

We also examined the continuous coherent state wave functions, obtained by VhLhS{z — z'), and found that these 
indeed seem to minimize the spatial dispersion. The CCS turn out not to be entirely homogeneous with respect to 
the guiding center position w but rather periodic with periods dictated by the lattice of broken translations xi x txi. 

We considered the effect of projecting a holomorphic derivative to the LLL. There we found that the projected 
derivative becomes a sum of translation operators T^llli^z ~ "D = '^lait^.. The projection could be used to sys- 
tematically transform derivatives to 2? but is computationally expensive, and does not necessarily respect modular 
invariance. 

We also proposed a map Vhcs that could be used instead of, or in conjunction with, 9^ — > 2? to map wave functions 
to the LLL. For simple hierarchical states, like those for = | where there are only particle condensates, the many 
body states already reside in the LLL and so Vlcs is redundant. For more complicated states, like — ^ where there 
are hole condensates, there naturally arise anti-holomorphic components of the trial wave functions. Here usage of 
'Plcs could be one way of mapping these functions to the lowest Landau level. 

Finally we considered an explicit example of states that have anti-holomorphic components, constructed trial wave 
functions for them in the coherent state basis and compared them numerically with the ground state of the coulomb 
potential. 
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Appendix A: Orthonormal basis states on arbitrary t torus 

We will here summarize some information on the basis functions that exist on the torus. These can be built from 
the cylinder functions under the constraints that 

t[Vn,s=e-'y^^^Vn,s (Al) 

*2f?n,s = rjn.s^k- (A2) 

As in the main text we parametrize our torus as Lx x L^t where L^t — La + iLy. Thus we think of the height of the 
parallelogram as Ly and the skewness as La- The flux relation is as usual La;Lj, = 2ttNs. The cylinder wave functions 
in the n:th Landau level have the form 

Xn,s - J— e-^y^-H^iy - y,)e-^(^-^=)'. (A3) 
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Under small magnetic translations they obey 



4x 



2 An,s 



e " Xn,s^k- 



to choose tts^t to be a p. 
that 



We see that to fulfill ( Al I we must construct rjn s — J2t '^s,tXn.s+tNs where j is a weight. To fulfill ( A2 1 we need 
lase. The linear combination needed has the form r/n^s = X^t 6*'^=+*"= Xn,s+tNs ■ We now have 



and want to choose Ps+tN^ such that Ps^tN^ ~ Vs+tN^'-^k + = Ps+tN^-k modulo 27r. We see that this is fulfilled 
if we choose /J^ = \vr^r such that 



(ys+ti»)(ws+tiA) 



The full wave functions in the LLL are thus 



(A4) 



which we can, up to a factor of e *2J's"Sj rewrite as 



e 2 



Appendix B: Cylinder functions on the torus 

Let us evaluate the overlap between the cylinder wave functions from ^ taken over the torus geometry. Normally 
one would not look at this since these functions do not form an orthonormal set on the torus. In our case we will be 
interested in this overlap, since when dz acts on periodized functions it will spoil the boundary conditions. To keep 
things general we will integrate over a torus resting with its lower left corner at a point z = a + t/3. The overlap, 
which is a dependent, is 

(Xm.|Xn.)torus='5../r'" (Bl) 

where 

IT " = ^ f ''^"^^ dy H„,{y)H„{y)e-y\ (B2) 

The integral has the property that X^t ^I+tAr ~ ^mn- With the knowledge of the overlap of cylinder wave functions 
on the torus we see that the rjn^s form an orthonormal set of states 

iVmr hns) = ^ 6^"'-+*"= {Xm,r+tN, IXn,.+iwJ torus = Ks^m^n- 

t,l 

From this follows that Vm — \ri„ir) {Vmr \ is a projector onto the mith Landau level and I — V„i is the identity 
operator. 

The rjn^s were given by rjn^s = X^r ^^^^''g*^^' "'^Xn.i-- Remembering that dxXn.s = —''•ysXn,s we see that we can write 
a dx derivative on rjns as 

t 

t 
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We also see that if we want to calculate the overlap of {r]m,r \dxiln,s) we end up calculating 

= —111 S S —iL p~'^hyr+lNs'^r+lNefpthVs+tNs'^3+tNs , , . ^ {-y , . , , \ 

— ''ys"m,n'Jr,s '-^y / \Xm,r+lNs \Xn,s+tNg 

Lt 



It is the part inside the parentheses that we will call G™'" = ysSm,n + Ly H^tN ■ t)y inspection that G™'" 

is periodic with G™'" = G^^^^ . 

Appendix C: Rewriting 

If we are looking at functions in the LLL we can simplify the expression of G^ = G°°. To facilitate further 
manipulations we will rewrite G^. The main difficulty comes from Ylit ^^s+tN which can be rewritten as 



t t Jct + tLy 



2 

where a = (5 — yg. This gives G^ as 



tL„+a ^ 

dy 

tLy—OL 



Gs=ys-^ + ^j: dye-y' (CI) 



which is numerically simpler since we are only integrating over patches of Gaussian functions. Reparametrizing 

2 



^" + d, such that 6 is the deviation from the symmetric integration region we have 



Gs=ys-^Y. dye-y . (C2) 
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